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GAMMA RADIATION FROM THE PROTON BOMBARDMENT 


ABSTRACT 


The nuclear gamma radiation following the non-resonant capture proton 
has been studied with protons energies from 800 kev. 2.1 Mev. and 
found consist three components. The main transition goes, via gamma ray 
energy that varies with proton bombarding energy, the state 
This state, 487 kev. above the ground state, radiates directly the ground 
state. addition there weaker direct radiative transition from the capture 
configuration the ground state. energies above 1.8 Mev. oxide 
targets bombarded with protons give rise radiation 873 kev. attributed 
the presence the isotope via the reaction the radiation 
corresponding the transition from the first excited state the ground 
state. 


INTRODUCTION 


Hitherto there has been report the detection nuclear gamma radi- 
ation accompanying the capture protons oxygen, although the produc- 
tion the positron active isotope the reaction has been 


Fi? 


observed several investigators (5, 10). Energy levels have been 
determined Laubenstein a/. (10) who measured the yield the elastically 
scattered protons from oxygen and the capture cross-section measuring the 
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activity the produced proton bombardment. assignment angular 
momentum the levels revealed has been made Laubenstein and 
Laubenstein (9). addition information has been obtained from study 
the neutron spectrum produced the reaction (1, 7), which estab- 
lished the existence excited state about 0.5 Mev. above the 
ground state. The energy level diagram for and for its mirror nucleus 
with presently accepted angular momenta given the well-known Energy 
Levels Light Nuclei, (2) and part this reproduced Fig. 

this work oxygen has been bombarded with protons energies ranging 
from 0.8 2.1 Mev., obtained from the U.B.C. electrostatic generator, and 
the non-resonant capture gamma radiation has been detected and measured. 
This has been found correspond the transitions marked and Fig. 
well the transition, from the first excited state the ground state. 


EXPERIMENTAL 


(a) Targets and Background 


Oxygen targets two sorts have been used, ice deposited gold backing 
and tungsten oxide formed heating thin tungsten sheet atmosphere 
pure oxygen. Both the pure gold and pure tungsten were cleaned electro- 
lytic etching followed careful washing, and the radiation 
from both these backings was found adequately small the region 
gamma ray energies main interest. typical background spectrum 
included Fig. was found that the spectrum and intensity from the 
tungsten and from the gold were nearly identical the gamma ray energy 
range 500 kev. Mev., and were little different from those the cosmic ray 
time-dependent background. The gold gave some extra radiation the 
Mev. region, about four times the tungsten background, while Coulomb 
excitation the low lying ‘rotational’ states gave rise considerably more 
background intensity and below 250 kev. for tungsten than for gold. 

(b) Detection 

The gamma radiation was detected means scintillation counter. Two 
were used, one consisting in. diam. in. high cylindrical block 
sodium iodide mounted Dumont 6292 multiplier and the other in. 
diam. in. high block mounted 6262 multiplier. The counter 
was carefully and heavily shielded reduce radiation from the electrostatic 
generator. Resolved proton beam currents were used, the beam 
being defined gold stop keep radiation minimum the target 
region. 

The pulse output from the counter after suitable amplification passed through 
cut amplifier and was displayed channel ‘kicksorter’. The kicksorter 
was set with precision, variable pulse generator and the energy scale for 
the counter was established frequent intervals using several standard sources 
for calibration. careful check the linearity the system showed 


better than over the energy range 0.51 2.62 Mev. 


WARREN PROTON BOMBARDMENT OXYGEN 


III. RESULTS 
(a) The Decay 

thick ice target was bombarded for several minutes and analysis made 
the gamma radiation emanating from the activity produced the ice. 
evidence for any gamma radiation, other than the annihilation radiation, was 
found. The background above 0.51 Mev. this case was very low and any 
nuclear gamma ray energy above 0.51 Mev. present the extent 
the annihilation radiation would, believed, easily have been detected. 

rather rough measurement was made the half life the activity pro- 
duced following the decay over several half lives using linear rate meter 
(4), the output which was fed Brown recorder. value 62.5 
was obtained accord with recent values (2) for the half life signi- 
ficant amount any longer lived activity was detected. 


(b) The Capture Gamma Radiations Variable Energy 

The gamma ray spectrum observed with protons energy from 0.8 1.7 
Mev. consisted three components represented 71, the energy 
level diagram Fig. typical spectrum, covering the gamma ray energy 
range from 0.6 2.4 Mev., obtained with fairly thin target tungsten oxide 
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Fic. counter pulse height distribution due radiation from tungsten 
oxide target bombarded with 1.355 Mev. energy protons, together with background due 
tungsten backing and time-dependent cosmic ray background. 
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The background radiation from clean tungsten target and the time- 
dependent cosmic ray background taken under similar conditions are also 
shown. clear that two gamma rays are present, the weaker one, 
corresponding the direct transition the ground state and the stronger one, 
corresponding transition level about 487 kev. above the ground 
state, presumably identical the even state reported Ajzenberg (1). 
The energy both these gamma rays varies with the proton bombarding 
energy the manner expected from the energetics the reaction. The 
numerical value the gamma ray energy for the ground state transition was 
found accord with 


within the experimental error +15 kev., including the uncertainty 
the exact target thickness. was established using the well-known 873 kev. 
and other resonances the reaction, and using the generating 
voltmeter interpolate this energy region. 

The yield both and increased very slowly over the energy range 
studied, accord with previous measurements the yield (10), and 
resonance characteristics whatever were noticed. The ratio the yield 
measured 90° the incident proton direction remained about 10/1 
over the energy range studied. Angular distribution measurements were made 
with difficulty owing the low intensity, and all that can stated 
that major departure from isotropic distribution occurred any the 
energies studied. the other hand the distribution the radiation was 
almost the pure pattern for from 1.0 1.9 Mev. For example 
1.90 Mev. data taken only from 90° and uncorrected for the solid 
angle subtended the crystal (0.026 sphere) fitted curve the form: 


The absolute cross-section was deduced from the counter efficiency estab- 
lished 1.3 Mev. using standardized Co® source and extrapolated using 
source and source, thus obtaining three points efficiency/ 
quantum energy curve between which the variation could interpolated 
calculation. The ice target thickness laid down the water vapor dispenser 
was established measuring the shift the resonance 1381 kev. 
occurring the reaction when ice was laid down thin 
calcium fluoride target. taking the stopping power ice this proton 
energy 6.2 e.v. cm.? per molecule, the number oxygen atoms 
present given target thickness could deduced. The absolute cross-section 
for the production radiation 1.35 Mev. was found 
6+3 


(c) The Transition from the 487 kev. Level 

The presence third component energy around 0.5 Mev. was made 
evident, spite the presence positron annihilation radiation, examin- 
ing the spectrum the 0.5 Mev. region over various times irradiation the 
target shown Fig. Thus the annihilation radiation activity, the end 
the first sec. run made with constant beam current, will only have 
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Scintillation counter spectrum 0.5 Mev. region from reaction taken 
(i) during initial sec. bombardment; (ii) after equilibrium activity had been established, 
beam current remaining constant; (iii) after beam was cut off, due activity; together 
with (iv) calibration, annihilation radiation peak from 
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Scintillation counter spectrum the 0.5 Mev. region from reaction 


taken during initial sec. bombardment, together with the spectrum resulting after the 
subtraction the annihilation radiation contribution due the positron activity the 
This latter contribution was the spectrum shown, which was obtained the sec. 
following min. bombardment. 
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reached 30% its final equilibrium activity achieved after some min. 
bombardment. Runs made during such initial short times indicated clearly 
the presence prompt nuclear gamma ray, energy just slightly below 
0.51 Mev., the annihilation peak. 

Measurements the positron activity, made over half life immediately 
following irradiation constant beam current lasting some nine ten 
minutes, enabled the equilibrium positron activity deduced, and hence 
the fraction the counts due positrons the peak obtained during the 
initial period run ascertained. Thus Fig. shown the pulse height 
spectrum the 0.5 Mev. region obtained during the initial sec. bom- 
bardment, together with the spectrum resulting after the subtraction the 
positron contribution. Likewise calculation the equilibrium annihilation 
radiation activity from the activity following irradiation showed that 
there were about 2.0 annihilation quanta for each 487 kev. quantum, compared 
the expected value 2.2 the basis branching ratio about 
10/1. This small discrepancy can accounted for the slightly higher coun- 
ter efficiency for 487 kev. than for 510 kev. radiation. 

Measurements made with thin target showed that the energy did not 
change with There was evidence for any specially large width this 
transition. The energy was below 0.51 Mev., the average value determined 
being 487 kev. This compared with the figure kev. 
given Ajzenberg from analysis the neutron groups from 
(1). would seem very likely that the level measured the same both cases, 
since the mirror nucleus has only one corresponding low lying state. 
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Fic. counter spectrum due radiation from oxygen bombarded with pro- 
tons 1.898 Mev. energy, background subtracted, showing presence 873 kev. component, 
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The angular distribution this prompt 487 kev. radiation was found 
isotropic. 

(d) The 873 kev. Radiation Which Appears for 1.8 Mev. 

Spectra taken proton energies around Mev. showed the presence 
fourth gamma ray, Fig. This gamma ray appeared with both tungsten 
oxide and ice targets and could not traced any impurity the target 
contamination the target from the vacuum system. 

Measurements showed that the energy this radiation did not change with 
and that was not detectable until was raised above 1.8 Mev. Its inten- 
sity rose sharply the proton energy was increased the maximum 2.1 
Mev. used this investigation. This radiation ceased soon the beam was 
turned off. The average value determined for the energy this gamma ray 
873 kev. Its angular distribution was found isotropic 1.9 
Mev. 

IV. DISCUSSION 


The absence any prominent nuclear gamma radiation from the decay 
accord with the beta ray data Wong (15) and indicates that 
only very few none the transitions via the 875 kev. level 
would expected from the ‘favored’ ‘allowed’ character the transition 
between the 5/2+ ground states and 

The gamma rays emitted during the capture proton fit the level 
scheme and energetically agree, within the limits error, with the binding 
energy the proton deduced from the threshold 1.836 0.003 
Mev. for the reaction measured Bonner and Butler (3). The 
experimental uncertainty the measured gamma ray energies considered 
about +10 kev. but addition there some uncertainty the exact 
value due uncertainty target thickness; the over-all error should not 
exceed +15 kev. 

The isotropic distribution the 490 kev. gamma radiation accord with 
the previous assignment for the first excited state (1, 2). 

However the fact that most the transitions are found cascade 
through the state 487 kev. shows clearly that the non-resonant 
gamma ray yield and capture cross-section for the production cannot 
attributed toa broad state this low energy, suggested Lauben- 
stein and Laubenstein (9). Two other explanations are possible. Firstly, the 
ordinary compound nucleus theory, successful for dealing with resonant 
reactions, applied, then the broad states 4.35 and 4.73 Mev. 
(10) above the ground level may invoked account for this yield. For ex- 
ample the 3/2— state 4.73 Mev. would give angular distribution the 
yield still quite measurable 800 kev., this rather far-fetched 
explanation and would still not account for the large cross-section for elastic 
scattering protons low energies (9). second alternative and more satis- 
fying explanation attribute the process direct radiative capture, with- 
out the formation proper compound nucleus. There are several cases 
such direct radiative capture reactions among the lightest nuclei, for example 
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the process, the and T(p, reactions 
(8, 11, 12, 13, 14), and the yp’) reaction (16), and may well that 
some other cases non-resonant gamma radiation arising from simple capture 
which the cross-section changes smoothly and slowly with energy are 
this type. 

even conceivable that the sin? pattern for the distribution may 
arise from the absence spin-orbit coupling, the case the D(p, 
reaction which the proton spin plays significant role. 

There remains the question the origin the radiation 873 kev. 
designated There little doubt that arises from oxygen, and the measured 
energy value very closely equal that the transition between the first 
excited state and the ground state. The isotropic angular distribution 
accord with this assignment since the first excited state even 
state. seems plausible attribute the radiation the reaction 
spite the low isotopic percentage present natural oxygen, 0.04%, 
since such process might well have moderately large cross-section. (The 
reaction prolific reaction, that its inverse will large, 
The rapidly increasing yield may attributed the level 7.6 Mev. 

might mentioned that the variable energy gamma rays together 
with the Mev. radiation, would seem provide rather convenient way 
establishing the relative counting efficiency scintillation counter for gamma 


radiation. 
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REFLECTION FROM WIRE GRID PARALLEL CONDUCTING 
PLANE! 


ABSTRACT 


solution outlined for the problem plane wave incident obliquely 
parallel-wire grid which backed plane conducting surface. The electric 
vector the incident wave taken parallel the grid wires. The equivalent 
transmission line problem pointed out. shown that, certain cases, 
resistive wire grid will absorb all the energy the incident wave. 


INTRODUCTION 


There have been numerous investigations the electromagnetic properties 
grating thin parallel cylindrical wires composed conductive material. 
Such gratings have been employed shield regions from spurious electro- 
magnetic fields act reflectors antenna systems simulate reflect- 
ing plane surface. The first quantitative study was made Lamb (6) 
1898 who considered plane wave incident normally grid. showed 
that the spacing between the parallel wires was small, the reflection and 
transmission the waves could varied changing the wire radius. 

1914 von Ignatowsky (5) made very exhaustive analysis the scatter- 
ing incident plane waves single, infinite metallic grids. His formulas 
have been rededuced, extended, and applied other authors (3, 12, 13) 
since that time. very illuminating treatment has been given MacFarlane 
(8) who indicated that single grid can represented impedance 
shunted across infinite transmission line whose characteristic impedance 
equal the intrinsic impedance free space. showed that this impedance, 
for any grid spacing, was proportional where the wire 
radius and correction factor which function the angle incidence 
and the spacing 

the purpose this paper consider the case reflection plane 
waves from wire grid which parallel conducting wall. this case 
can anticipated that the grid may still represented shunt impedance 
whose value, however, will depend certain extent the distance 
the conducting plane. will shown this paper that the shunt impedance 
incidence, grid spacing wavelength ratio, and the distance the grid the 
conducting plane. The resultant reflection coefficient can then found directly 
representing the conducting wall short circuit across the equivalent 
transmission line mentioned above. 

The possibility modifying the properties the grid inserting small 
resistors the wires also discussed. appears that under certain conditions 
the reflection coefficient can made vanish the distance the conducting 
plane behind the grid properly chosen. 


received May 14, 1954. 
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GENERAL THEORY 


With respect Cartesian coordinate system the wire grid contained 
The grid composed array wires parallel the axis and spaced 
distance between centers. The wires are taken circular cross section 
and are assumed initially infinite conductivity. The insulating space 
surrounding the wire grid homogeneous and has dielectric constant 
and permeability 

plane wave with the electric field magnitude parallel the axis 
impinges the grid with angle incidence indicated Fig. la. The 
primary incident field given 
where the time factor has been omitted and 


where the wavelength. The currents induced the wires will all 
equal magnitude but will have progressive change phase between 


e 
FIGURE FIGURE 


Fic. The parallel wire grid and the conducting plane. 
The equivalent problem employing the image principle. 


can arbitrarily zero. The incident wave will also reflected from the 
conducting plane calculate the magnitude the current the 
wire and the resultant reflected fields convenient apply the theorem 
images. 

The conducting plane and the wire grid are replaced two identical wire 
grids contained the planes shown Fig. satisfy the 
boundary condition that the electric field should vanish the plane 
reflected wave has the form 


e e 
e e 
H, H, © e 
e e e 
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REFLECTION FROM WIRE GRID 


and addition the currents the image grid must equal and 

The problem now find the field point P(x, result scattering 
the two plane waves from the double grid system. assumed that the 
wires are small diameter compared with the spacing and with the distance 
distance from infinite wire carrying current well known (3) and 
given 


where the Hankel function the second type and chosen give 
rise outgoing waves infinity for time factor exp(iwt). The field 


which can transformed more convergent summation (see Appendix), 
follows: 


where Correspondingly, the field due the currents flowing 

The total field then 


current now found from the boundary condition that the field must 
also zero the surface the wires the grid. The point taken 


where 
with 
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where h/X, and has been assumed that >a. The current 
incident field 

(no /2 cos — + 


This result along with equations [3] [6] constitutes the solution the 
problem. The restriction infinite conductivity the wires composing the 
grid can removed the boundary condition the wires replaced 
the condition 


where the internal impedance the wire. can assumed that the 
field uniform around the wire since and <h, and hence can 
calculated known methods (10) and given 


[12] 


permeability, conductivity, and dielectric constant the wire material. 
and are modified Bessel functions order zero and unity. For metallic 


wires the displacement currents are negligible since even microwave 
frequencies. addition, the frequency usually sufficiently high that 


Employing this result follows that the defining equation for should read 


which reduces equation [8] for the case 


THE DISTANT FIELD 


The case particular interest that which the field point sufficiently 
far removed from the grid that inspection equations [3] 
and [4] for the field and evident that only the term 
significant for d/\ The higher values correspond 
evanescent waves which are highly damped the positive direction. For 
larger values additional undamped waves can scattered from the 
grid. The discussion will limited here, however, the smaller grid spacings 
satisfying the above inequality. The field can then written follows: 


cos 


WAIT; REFLECTION FROM WIRE GRID 


Employing equation [10] this written 


where 


[16] 


Evidently the first term [16] the reflection coefficient plane conducting 
surface whereas the second term represents the effect imposed the presence 
the wire grid. 

Fig. illustrates electrical circuit which may taken the analogue 
the wire grid and conducting surface. 


Fic. The analogous electrical transmission line. 


The space the right the plane represented transmission 
line whose characteristic impedance and whose propagation constant 
distance from the shorted end impedance shunted across the line. 
The voltage considered due the superposition the two travelling 
waves, that 


and the current therefore 


can verified readily from transmission line theory that given 
equation [16] 


and are then identified with the electric field and the magnetic field 
component respectively. 

The factor can now termed the shunt impedance the grid and 
given equation [13]. can seen that for and the factor 
negligible compared with and therefore 


which identical with the expression for the impedance single resistive 
wire grid free space for the case and normal incidence given 
Gans (2). 
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apparent therefore that the grid possesses shunt impedance which 
mainly function and wire radius and which depends lesser 
extent the angle incidence and the presence the conducting surface. 
Values the correction factor are plotted Fig. (d) function 


Fic. correction factor function the distance the grid the conducting 
plane for specified angle incidence. 


APPLICATION 


possible application this analysis the design simulated anti- 
reflective coating, for letting 
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[20] =iX, +R, 


where and are the shunt reactance and shunt resistance the grid, 
can seen easily from equation [16] that then 


and 


where equation [22] immediately implies that the wires must possess finite 
conductivity otherwise modified absorb the energy the incident 
wave. 
The equation for the reactance, for the case 1/(1 sin then 
written follows: 
(no/2 cos 


apparent that the sine must negative, therefore lies the range 

For conducting grid fabricated from metal wires the actual resistance 
will usually much smaller than specified equation [22]. should 
relatively simple matter, however, modify inserting small carbon 
resistors ohms regular intervals along each the wires. For 
and can seen that 


where the second term, which the effective resistance the wire, would 


small compared with the first term. The required value for vanishing 
reflection coefficient then given 


Equations [23] and [25] can employed calculate the required value 
and for zero reflection having specified values and When 
the grid spacing becomes small, however, such that the design equa- 
tions become much simpler. this case cos slightly greater than odd 
number quarter wavelengths slightly less than integral number 
half wavelengths. For example 


for 1,2,3..., then and are small compared with since 
small quantity compared with unity this order approxima- 
tion follows from equations [22] and [25] that 
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2 = = 
and 
for the conditions equation [26] and 


for the conditions equation [27]. 

experimental study wire grids backed conducting and dielectric 
media being carried out this laboratory. The results will reported 
later date. 


APPENDIX 


Summation Formula 
derivation given here for the relations expressed equations [3] and 
[4]. The starting point Poisson’s formula 


where 
+oo 
[32] H(m) h(x) 


This formula valid h(x) integrable the interval con- 
tinuous, and steadily decreasing zero and H(m) the Fourier 
transform h(x). Letting 


where use has been made formula Campbell and Foster (1) (No. 917 
now given 


+00 
m=—co k)° a’ | 


When set equal zero the above equation reduces formula derived 
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Smith, Infeld, and Chien (11). When also allowed approach zero 
follows that 


where log ..., Which Euler’s number. The special case 
equation [36] where was given originally Ignatowsky (5). 
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THE PHOTONEUTRON AND PHOTOPROTON CROSS SECTIONS 
SILICON AND MAGNESIUM! 


ABSTRACT 


The (y, and cross sections were measured function photon 
energy number isotopes. They were found have the resonance shapes 
characteristic photonuclear reaction for which the pertinent data are sum- 
marized the table. 


Maximum Integrated 
Position cross cross 
Isotope Reaction Threshold of peak section section 
(Mev.) (Mev.) (mbarns) (Mev-barn) 
Mg® (y, 7.2 ~13 ~l4 ~0.05 
(y Pp) 12.1 20.5 27 0.12 
Meg* (ym) 11.1 ~17 ~30 «0.09 
(yp) 14.3 22 25 0.07 
Si? (y, n) 8.5 and 10.6 15 23 0.096 
+Si” (y, Pp) 12.3 and 12.9 2 2 0.23 


view the fact that each isotope the threshold for the (y, reaction 
higher than that for the (y, reaction surprising find the maximum 
(y, cross section least equal the maximum cross section, 
while peaking about Mev. higher energy. These results are discussed some 
detail. 


INTRODUCTION 

The applicability statistical nucleus theory (32) photo- 
nuclear reactions the energy region Mev. still much question. 
Considerable evidence may found the literature for its apparent failure. 
For instance, the energy spectrum protons emitted photodisintegration 
experiments seems indicate high energy tail (5, 31) not expected from the 
statistical theory. Measurements the angular distribution the emitted 
protons and neutrons seem indicate anisotropy some cases (4, 12), 
fact which not easily reconciled with statistical theory. Again, experiments 
Ziirich (11) with y-rays from the reaction energy 14.8 and 
17.6 Mev. indicate greater ratio proton neutron reactions than pre- 
dicted the statistical theory. all these cases however the deviations from 
the statistical theory involved relatively small fraction the total number 
reactions occurring. For example, the number protons the high energy 
tail (5) the proton spectrum relatively small fraction the protons 
involved. 

explain these discrepancies Courant (3) and Marquez (20) proposed that, 
addition the usual (8) type interaction which always 

received May 31, 1954. 
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results compound nucleus, some cases photon absorbed one 
proton the nucleus, which then emitted without the formation inter- 
mediate compound nucleus. According Courant the protons emitted the 
direct process are usually the statistical protons. They should 
most easily observed nuclides whose threshold not much lower than 
the (y, threshold. This precisely the case the elements examined 
Hirzel and (11), and relatively small contribution from the direct 
process sufficient explain their results Courant has shown. must 
pointed out however that the measurements Hirzel and refer 
(y, and processes different isotopes the same element. has been 
found that strong fluctuations may occur the photo cross sections such 
cases and their comparison must interpreted with some care. Similarly 
small contribution from the direct process will serve explain the high energy 
proton tail observed Diven and Almy (5). 

would desirable compare and cross sections the same 
isotope. Outside the measurements Halpern and Mann (9) and Toms 
and Stephens (31), very few cross sections are available for such com- 
parison. These workers have published the (y, cross sections for 
these elements have been measured our laboratory (21, that inter- 
comparison possible these cases. turned out that the measured ratio 
Was quite satisfactory when compared with that expected from 
statistical theory, except possibly for the case analysis the 
thresholds for these elements brought light the fact that all cases examined 
the threshold was lower than the (y, threshold least Mev. 
Some these thresholds are given Table 


TABLE 
THRESHOLDS FOR AND REACTIONS THE ISOTOPES INDICATED 


Element (y, Ref. (y, Ref. 
15.95 (19) 18.7 (19) 
8.26 (18) 12.75 (28) 
7.64 (13) 10.25 (28) 
7.5 (9) (13) 
8.8 (9) 11.6 (13) 
6.3 (9) 8.7 (28) 


This particular sequence thresholds favors the emission protons the 
excited nuclei. Protons resulting from the direct process would not this case 
easily detected. immediately seen that examination some iso- 

which the threshold higher than the (y, threshold would 
great interest. Fortunately, the low percentage isotopes magnesium 
and silicon present two such cases, can seen from Table II. For this 
reason these cross sections have been measured and the results are presented 
this paper. 
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TABLE 
PHOTONUCLEAR THRESHOLDS THE MAGNESIUM AND SILICON ISOTOPES (18) (MEv.) 


Isotope Abundance, (y, (y, 
78.8 16.6 
10.1 
92.2 11.6 16.9 
Si?9 4.7 12.3 8.5 
12.9 10.6 


EXPERIMENTAL PROCEDURE AND RESULTS 


(a) Silicon 

The reactions and lead 2.3 min. half-life, and 
6.6 min. half-life. One-gram samples silicon powder 99.87% pure were irradi- 
ated cadmium thimble inside lucite block. After irradiation the decay 
curves the silicon powder were measured and resolved into 2.3 and 6.6 min. 
components. other activity was observed. 

Absolute determination the activities was carried out Mev. using the 
reaction standard. Samples thick enough give cosine 
distribution the emitted rays were used and the procedure outlined 
Baker and Katz (1) was followed. The decay scheme used for was that 
Motz and Alburger (22) and for that Seidlitz, Bleuler, and Tendam 
(26). The uncertainty the branching ratio the disintegration scheme 
will result uncertainty about this activation curve. The 
absolute values the cross sections are thought accurate within 20%, 


Activations 


Peak Energy Mev Peak Bremsstrahlung Energy Mev 


Fic. Activation curves resulting from Fic. Neutron yield 
(y, reactions and silicon. 
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RESIDUAL 


Mev 


Fic. cross section natural silicon. The dashed curve represents the 
(y, cross section corrected for isotopic abundance, measured Summers-Gill 


al. (29). 
| 
Residual 
Activity 


Mev 


Fic. cross sections and The (y, curve has been normalized per mole 
these isotopes mixed according their isotopic abundance. The (y, curves are per mole 
the respective isotopes. 
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which the figure usually quoted our measurements and which has been 
substantiated measurements other laboratories. The relative cross 
sections are probably known with somewhat greater precision. The positions 
peak cross section are established within Mev. 

The activation curves are shown Fig. they were solved for the cross 
sections using the difference (15). The results are shown 
Fig. 

The (y, reactions cannot measured residual activity since the end 
products the reactions are stable. For these measurements therefore used 
method which detected the emitted neutrons. This method has been pre- 
viously described (21). has the inconvenience giving the total cross section 
for natural silicon. Advantage was taken, however, the fact that the threshold 
for the (y, reaction the main isotope, very high (29) (16.9 Mev.). 
Thresholds for the other isotopes are much lower the neutron yield curve 
exhibits the peculiar shape shown Fig. Normalization this curve was 
carried out against copper number energies. The cross section obtained 
shown Fig. The cross section obtained Summers-Gill al. (29) for 
the reaction properly corrected for isotopic abundance shown 
the same figure dashed curve. There excellent agreement between both 
curves above the threshold for the (y, reaction hence the little 
and separation the components here possible. 


(b) Magnesium 


measure the cross sections magnesium, thin disks this material 
were irradiated the betatron beam various energies and counted after the 
sec. activity resulting from the (y, reaction had died out. The 
decay curve obtained had sec. component from the 
reaction and 14.9 hr. component from the reaction. The 
activation curves were normalized the known yield Mev. 
after making the necessary corrections for self-absorption (1) and decay scheme 
(23). 

These normalized activation curves are shown Fig. They were analyzed 
the photon difference method and yielded the cross section curves shown 
Fig. Comparison these cross sections previously published values will 

the case silicon, the (y, cross sections were measured detecting 
the neutrons resulting from these reactions. The neutron yield curve, stand- 
ardized against the neutron yield from the reaction copper according 
our usual procedure (21), shown Fig. Analyzing this curve the 
photon difference method resulted the cross section curve shown Fig. 
This cross section course the sum the cross sections the three 
isotopes; however the large peak above 16.6 Mev. must due the 78.8% 
isotope which has its threshold this energy. 

The cross section for this isotope had been previously measured Katz 
and Cameron (16) detecting the sec. activity resulting from the 
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Peak Bremsstrahlung Energy Mev 


Fic. Activation curves resulting from (y, reactions and Note the scale 
the right-hand side for the curve. 


activations mole 


Peak Bremsstrahlung Energy Mev 


Fic. Neutron yield from natural magnesium. 
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26 
activi 


6 8 10 12 14 ie. 18 ao 
Mev 


Fic. cross section natural magnesium (solid line). The dashed curve for 
was obtained from the (y, activation curve Fig. corrected the isotopic abundance 
the natural element. The difference between these curves represents the cross sec- 
tions the 10.1% isotope and 11.1% isotope 


© Katz and Cameron 
+ Remeasured points 


activations mole 


12 14 16 20 22 24 26 28 30 
Peak Bremsstrahlung Energy- Mev 


Fic. curve representing the (7, reaction normalized according 
the procedure outlined the text. The remeasured points were taken during the course the 
present experiments check the older data Katz and Cameron (16). 


reaction. Since rather difficult normalize properly the 
sec. activity, their result 9.8 mbarns for the peak cross section somewhat 
inconsistent with the curve Fig. order obtain better value 
normalized their activation curve value consistent with the neutron yield 
data Mev. from the natural element. This assumes that good approxi- 
mation far above threshold the separate isotopes have equal neutron yields 
per mole per roentgen. Since forms about 80% the natural element, any 
our present experiment the neutron yield Mev. was found 
1.75 neutrons/mole/100 This excellent agreement with 1.64 107 
found Price and Kerst (25) and 1.72 107 Montalbetti, Katz, and Gol- 
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demberg (21). thus normalized the (y, curve the average value 
1.70 activations/mole/100 This activation curve shown Fig. 
check repeated some the points which are also shown this figure. 
The cross section obtained from this activation curve normalized for isotopic 
abundance has been superposed Fig. 

the systematics (y, reactions previously found (21) apply these 
reactions then one would expect the peak cross section for each the isotopes 
lie approximately Mev. above threshold, except perhaps the case 
which because its very high threshold might expected behave 
more like which the peak cross section only about Mev. above 
threshold. This prediction fact borne out the cross section curve Fig. 
resolved indicated the dashed lines. The initial bump the curve 
presumed predominantly due the cross section, while the 
difference between the cross section and that obtained for 
residual activity peaks about Mev. above the threshold 
and presumed the cross section for this reaction. 

Though the method used resolve these cross sections has been tested 
experimentally and found quite accurate (6), conservative would 
suggest that the peak cross sections have been established absolute value 
within factor and the positions the peaks within Mev. These 
two resolved cross sections corrected for isotopic abundance have been replotted 
Fig. where they are compared the cross sections for the same 
isotopes. 


Activity 
M 


Katz Cameron 


6 8 10 12 14 16 18 20 22 24 26 


Fic. Photo cross section and The cross section published 
Katz and Cameron (16) shown the dashed curve. 
DISCUSSION 


(a) Silicon 
The characteristics the reactions silicon measured these experi- 
ments, corrected for isotopic abundance, are summarized Table 
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TABLE III 
CHARACTERISTICS THE PHOTONUCLEAR REACTIONS THE SILICON ISOTOPES 


Width Maximum Integrated 
Position half cross cross section 
Reaction peak maximum section Mev. 
(Mev.) (Mev.) (mbarns) (Mev-barn) 
20.5 3.5 0.075 


Since not possible resolve the cross sections the and 
isotopes have combined the cross sections these isotopes according 
their isotopic abundance. This averaged cross section shown 
the last line the table and compared the average cross 
section shown the second line. 

number surprising facts are immediately evident from such compari- 
son: (a) the peak cross sections the two curves are Mev. apart, the 
(y, curves are much wider, and (c) the peak cross sections and integrated 
cross sections for the reaction are considerably larger than for the (y, 
reaction. This last fact particularly surprising since the (y, threshold 
these isotopes considerably higher than the (y, thresuold. For better 
comparison, the two curves are drawn the same diagram Fig. 
immediately evident from this diagram that the two cross sections are about 
equal Mev. whereas statistical theory would predict 0.05. 
The divergence between theory and experiment even greater higher energy. 

view the above conclusions our results have been examined with great 
care for possible errors. The results are deemed quite reliable since 
they were taken relative the reaction, which was measured 
Summers-Gill, Haslam, and Katz (29) using residual activity, and indepen- 
dently the present experiment using the neutron detection method. The two 
measurements are excellent agreement and good agreement with the (y, 
cross section systematics established Montalbetti, Katz, and Goldemberg 
(21). 

establish the accuracy the (y, measurements particular signi- 
ficance point out the fact that the sum the integrated cross section the 
two isotopes for both reactions 0.33 Mev-barn whereas 
the theory (23) would require minimum integrated cross 
section (minimum provided their formula) 0.44 Mev-barn. This would 
indicate that our results are certainly not too large, and since the (y, con- 
tribution has been quite well established, the integrated cross section 
certainly not too large. Furthermore, the value 0.33 Mev-barn for the sum 
both reactions falls smoothly the curve defined the measurements 
Montalbetti, Katz, and Goldemberg (21) carbon, aluminum, cobalt, and 
niobium. The only values the literature the (y, cross sections the 
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silicon isotopes the energies involved are those Hirzel and (11). 
Unfortunately their results are given relative the (y, cross section 
There wide discrepancy between their measurements this latter isotope 
and those made our laboratory, that comparison for the silicon isotopes 
would quite meaningless. 

Perlman and Friedlander (24) using and 100 Mev. bremsstrahlung have 
measured the ratio the proton yield from the neutron yield from 
They find this ratio 6.6 Mev. and 5.8 100 Mev. these energies 
the ratios the yields photonuclear reaction may expected equal 
the ratios the integrated cross section. Johns al. (14) report the inte- 
grated cross section the peak 0.015 Mev-barns. assume 
that this cross section follows the usual shape and symmetrical about the 
peak, the total integrated cross section 0.03 Mev-barns. Using 
this value and Perlman and ratio estimate the integrated 
cross section ~0.18 Mev-barns, excellent agreement with the 
value 0.19 Mev-barns reported here. 

Sheline (27) has measured the ratio the yield relative the 
(y, yield 35-48 Mev. and found this 2.7 0.3. Following the 
procedure the previous paragraph and using 0.05 Mev-barns for the inte- 
grated cross section (10), the integrated cross section 
estimated ~0.13 Mev-barns. The agreement here not quite good but 
satisfactory considering the errors involved. Sheline also measured the ratio 
+0.2. This agrees favorably with the ratio the integrated cross sections for 
these reactions which find 1.4 0.2. 


(b) Magnesium 
The characteristics the reactions magnesium measured these 
experiments, corrected for isotopic abundance, are summarized Table IV. 


TABLE 


Width Maximum Integrated 
Position half cross cross section 
Reaction peak maximum section (Mev-barn) 
(Mev.) (Mev.) (mbarns) 
19.5 5.8 14.2 0.093 
20.5 3.5 0.12 


The values the (y, reactions and are indicated approxi- 
mate according our previous discussion and may taken indicate that 
the peak positions have been determined within Mev. and the cross 
sections better than factor Other values are estimated within 
our usual accuracy +20% for cross section and Mev. for peak position. 
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spite these wide limits immediately evident that the peaks the 
(y, and cross sections the same isotope are widely separated (an 
average Mev.) the case silicon, and that there energy region 
for each isotope where the ratio greater than one and increasing 
function energy. 

view these conclusions will examine whether the data Table 
are consistent with our knowledge photonuclear reactions and previously 
published measurements. The element nearest mass number magnesium 
whose photonuclear properties have been measured aluminum. Halpern and 
Mann (9) found the integrated cross section this element 0.12 
Mev-barn and measurements our laboratory gave the integrated 
cross section (21) 0.045 Mev-barn. Therefore the total integrated cross sec- 
tion for aluminum 0.165 Mev-barn. This excellent agreement with 
0.17 Mev-barn for and 0.16 Mev-barn for obtained from Table IV. 
Since integrated cross section expected smooth function mass 
number, the agreement the above values gratifying. 

The (y, and cross sections cannot compared separately those 
neighboring nuclei, since the ratio strongly dependent their 
relative thresholds and the neighboring nuclei have the thresholds lower 
than the thresholds. The (y, cross sections have already been discussed 
and are presumed have been established better than factor Measure- 
ments Katz and Cameron (16) magnesium gave the cross section 
shown the dashed curve Fig. Except for initial flat portion 
the newer curve both cross sections are excellent agreement within 
experimental accuracy. For Katz and Cameron found peak cross section 
14.8 mbarns 21.7 Mev. and Toms and Stephens (30) measured the same 
cross section using nuclear plates and found have peak value 13.7 
mbarns 20.5 Mev. with uncertainty the order 50%. Our present 
measurements give mbarns and 20.5 Mev. for these values. The positions 
the peaks are good agreement all three measurements, but the newer cross 
section considerably higher. The reason for this discrepancy present not 
known; however examination Fig. shows that even the 
cross section were reduced factor and the cross section 
were increased factor would still not possible reconcile our 
results with the usual assumption competition between neutron and proton 
emission from excited nucleus. 


QUALITATIVE ANALYSIS 
Our results may explained tentatively the assumption that the process 
proposed Courant plays dominant role photonuclear interactions. The 
following arguments will show that this assumption not necessarily contrary 
any previous experimental results and will satisfactorily account for our 
observations. shall discuss the usual (8) and Courant (3) 
interactions separately. 


Interaction 
this case the nucleus whole interacts with the incident gamma-ray 
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quantum forming intermediate compound system. The excited nucleus 
then decays manner suggested Weisskopf’s statistical nucleus model. 
Courant Interaction 

The incident gamma-ray quantum interacts with single proton the nuc- 
leus, this proton then either escapes giving rise the direct photoprocess 
loses its energy collision within the nucleus. this latter effect 
occurs compound nucleus formed which indistinguishable from case 
observe the direct photoprotons important suppress process 
and favor (B,); furthermore important. that the proton/neutron ratio 
expected from small that the observed protons may properly iden- 
tified coming from 

The proton/neutron ratio from small when the threshold higher 
than the (y, threshold, which precisely the case the isotopes investigated. 
Process reduced importance when the coulomb barrier favoring 
the escape the direct protons. equal importance that the nucleus have 
large surface volume ratio that the proton escape enhanced; further- 
more, necessary that the nuclear diameter small that the initial 
motion the direct proton into the interior the nucleus, its mean free 
path may commensurate with the nuclear diameter. All these require- 
ments are admirably satisfied the magnesium and silicon nuclei. Previous 
measurements have not been conclusive since relatively large number 
protons were expected from and particular since process was important. 

Recent measurements Butler and Almy (2) the (y, and 
cross sections and will now examined the light the above 
discussion. For they find that the (y, threshold 2.4 Mev. higher than 
the threshold; however, since this relatively large nucleus, process 
important that the cross section does not exceed the cross 
section though their ratio about 100 times larger than expected from 
statistical model calculation. For the (y, threshold higher than the 
threshold 5.1 Mev. this case the number protons from processes 
and expected quite large though again small but definite depar- 
ture from the statistical model found. 

Calculations the shape expected for the direct process (y, cross 
section indicate that should have resonance shape. This results from 
initial portion which rises slowly because barrier effect and final decreasing 
portion falling approximately because the decreasing strength 
interaction between the photon electromagnetic field and the proton. very 
rough computation places the peak about Mev. from threshold the 
magnesium and silicon isotopes; experimentally the peak found about 
Mev. from the threshold. Marquez (20) has suggested that the peak should 
occur about twice the threshold energy. 

must again emphasize that the above discussion highly speculative 
and only presented this time indicate how our experimental results 
may fitted into coherent picture. 


The authors wish thank the National Research Council Canada for 
financial assistance. 
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ABSTRACT 


The (y, and (y, cross sections have been measured function 
photon energy from threshold Mev. The (y, cross section was found 
have peak value mbarns Mev. and integrated cross section 
0.35 Mev-barn. These values for the cross section were found 120 
mbarns >25 Mev. and 0.54 Mev-barn respectively. spite the higher 
(y, threshold, was found that the two cross sections were equal 19.5 Mev. 
and above this energy the (y, cross section became larger. This interesting 
result discussed. 


INTRODUCTION 

Recent work our laboratory, presented the accompanying paper 
(hereafter referred has again emphasized the importance studying 
the and (y, cross sections the same isotope for better under- 
standing the photonuclear process. Since the isotopes examined (I) are all 
low abundance their natural element, was felt that emphasize the 
conclusion reached would desirable study isotope high abundance. 

examination the light weight isotopes revealed that was the only 
one with the desired properties. has abundance 99.6% natural 
argon and its threshold (12.44 Mev.) higher than the (y, threshold 
(10.25 Mev.). Furthermore, the (y, reaction results the 55™ activity 
which may conveniently measured methods previously developed 
(6), and the (y, cross section can measured detecting the neutrons 
emitted the reaction using apparatus described Montalbetti, Katz, and 
Goldemberg (12). 


EXPERIMENTAL PROCEDURE AND RESULTS 


not possible detect this reaction residual activity measurements 
because the long half life However, using liquid argon samples, was 
possible detect the neutrons emitted the reaction with the aid the 
apparatus previously described (12). 

The apparatus used liquefy the argon gas shown Fig. Lique- 
faction resulted when the gas was cooled with the aid commercially available 
liquid nitrogen. This method very suitable for producing liquid argon the 
small quantities required. order prevent excessive pressure from building 
the container when removed from the Dewar flask, necessary 
open the valve the line the atmosphere just before doing so. Closing 
this valve and opening the valve attached the stainless steel container forces 


‘the liquid out. Experiments were performed with the liquid argon specially 


constructed thermos flask shown Fig. 
The X-ray beam was collimated that all passed through the liquid 


1Manuscript received May 31, 1954. 
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pressure gauge 
stainless steel 
container 


argon 
cylinder 


35mm. glass tubing 


beam 


seal liquid 
silvered argon 


Fic. (a) Method used for the argon gas. 
(b) Thermos flask for holding the liquid argon during irradiation. 


argon. Measurements could carried for about eight minutes before 
evaporation caused the top surface the liquid drop below the top edge 
the X-ray beam. During this time was possible measure least three 
points the activation curve (Fig. 2). One precaution had observed 


activations mole 


peak Bremsstrahlung 


Fic. Neutron yield per mole argon per 100 resulting from the (y, reactions, and 
the number activations per mole per 100 resulting from the (y, 


these measurements; since the boiling point argon below that oxygen, 
condensation oxygen from the atmosphere gradually replaces the argon 
the sample. keep this effect minimum replaced all the liquid the 
sample container after each set three measurements. 


1/16" monel tubing 
flask 
nitrogen 
valve 
(b) 
4 J 
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dry run with the empty container place was made determine the 
neutron background. This was subtracted from the total neutron yield give 
the curve shown Fig. Each point this curve actually the mean 
four separate runs. The ordinate given neutron yield per mole 
per 100 roentgens irradiation the indicated peak bremsstrahlung energy. 
This normalization was made comparison with the neutron yield from 
given copper sample 18, 20, and Mev. The yield from copper had pre- 
viously been well established (8) and confirmed independent measurements 
other laboratories (4, 10). The cross section was computed from Fig. 
the photon difference method (8) and shown Fig. 


10 12 14 16 18 22 @& 
Mev. 


Fic. Photoneutron and photoproton cross sections function photon energy. 


The (y, Reaction 


The yield from the (y, reaction was measured detecting the resulting 
activity. Argon gas pressure about two atmospheres was 
irradiated iron cylinder cm. long and cm. diameter with about 
1500 roentgens. The formed the (y, reaction chemically active and 
adheres the cylinder walls. After irradiation the argon gas was allowed 
escape and the chamber was flushed with 100 ml. 0.006 hydrochloric 
acid. The hydrochloric acid solution was then heated almost boiling and 
ml. 0.2 silver nitrate added, along with few grains solid detergent 
aid coagulation the silver chloride. 

The precipitate was centrifuged out, then suspended ethanol and filtered, 
the resulting sample being about thick. The activity was 
measured standard counting arrangement. Decay curves were plotted for 
several samples and these samples were found decay with pure 55™ half 
life. 

was assumed that all the formed the irradiation was recovered 
the above procedure. prove that chlorine was carried out the gas 
released after irradiation, the escaping argon one the experiments was 
bubbled through dilute hydrochloric acid. measurable activity could 
detected this solution. The exchange between chlorine atoms the washing 
solution and the radioactive nuclei adhering the cylinder walls after 
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irradiation is, course, very efficient because the numerically larger number 
chlorine atoms the washing solution. 

obtain the absolute yield from the (y, reaction the activity each 
sample was subject the usual correction for self and external absorption (1) 
well for back scattering and counting geometry. The mass each sample 
was calculated from the volume the irradiating chamber and the pressure 
and temperature the enclosed gas. Since the gas container subtended 
fairly large angle the X-ray target during irradiation, care was taken 
properly average the angular dependence the beam over its surface. this 
manner each irradiation was normalized the same dose rate. The yield 
radioactive nuclei per mole per 100 roentgens irradiation shown 
Fig. Each point this curve the average four separate experiments. 
The values below Mev. were reproducible within experimental accuracy, 
estimated better than 2%. Above this energy there was considerably more 
scattering, individual values varying much 10%; the average, never- 
theless, forming smooth curve. 

The cross section obtained from the activation curve shown Fig. 
Starting from threshold exhibits the slowly rising portion characteristic 
charged particle reactions until the excitation energy sufficient overcome 
barrier effects, then rises very steeply. Unfortunately, the maximum operat- 
ing energy our betatron too low reach the peak this cross section. 
Though would desirable also have the high energy end the cross 
section, the portion obtained sufficient for the present analysis. 


DISCUSSION 


The characteristics the cross sections are summarized Table Previous 
measurements (12) various elements Mev. indicate expected 


TABLE 
CHARACTERISTICS PHOTONUCLEAR REACTIONS 


(Mev.) (barns) Mev. 
(Mev-barn) 


0.038 


0.35 
>25 >0.12 0.54 0.8X 


measurements Ferguson, Halpern, Nathans, and Yergin, Bull. Am. Phys. Soc. 
29: No. 1954. 
23.5 Mev. 


neutron yield 10° neutrons per mole per 100 20. similar smooth 
interpolation previous results (12) would expect for the reaction 
peak cross section ~0.04 barn and integrated cross section ~0.3 
Mev-barn. These values are good agreement with those given Table 
Our results are also excellent agreement with the recent measurements 
Dr. Halpern Pennsylvania shown Table Proton yields from the 
reaction have been measured few elements only (5). The reported yields 
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protons per mole per 100 23.5 Mev. are 0.4 for aluminum and 
1.2 108 for nickel. Our yield 0.8 argon the same energy 
excellent agreement. The sum the two integrated cross sections 0.89 
This equal that obtained from the equation Levinger and 
Bethe (11) with though this region the periodic table the exper- 
imental values are usually about two thirds the theoretical value.? For 
greater ease analysis, the cross sections and their ratio are listed Table 
1-Mev. intervals. plot against (hv)! yields straight line, though 
the reason for this interesting behavior not present understood. 


TABLE 
energy (mbarns) (mbarns) 
2.8 
0.11 
3.4 0.17 
6.8 0.26 
0.60 
0.89 
1.26 
1.8 
3.6 
107 


examination Table confirms the conclusions reached (I) that 
difficult reconcile these data with statistical theory; above Mev. the 
(y, cross section increasing function energy while the (y, cross 
section decreasing function. important note that above 19.5 Mev. 
the ratio becomes greater than one that are not dealing with situation 
investigated number previous workers (2, 14). They were working 
with heavy nuclei and found that the ratio was greater than expected 
from statistical theory, though these nuclei Thus Toms and 
Stephens (14) have recently measured the proton yields from foils indium, 
barium, cerium, and bismuth when irradiated with Mev. bremsstrahlung 
X-rays. They found these yields ~1/100 the neutron yields previously 
measured Price and Kerst (13). The ratio predicted theoretical calcula- 
tions based evaporation from statistical nucleus was found them 
from ~1/3 ~1/100 that observed. explain these results Toms and 
Stephens analyzed the energy distribution the observed protons and con- 

that somewhat over 50% all those emitted are the result direct 
photonuclear effect. The over-all photonuclear process however not greatly 
altered these considerations, since mentioned above only ~1% the 
emitted particles are protons. 


for example Fig. Reference (12). 
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have pointed out that the evaporation model and direct photoeffect 
model for proton emission tend become indistinguishable for heavy nuclei 
with large volume surface ratio and high barrier for proton emission. 
Thus difficult separate these two processes the heavy nuclei used 
Toms and Stephens. Argon the other hand light nucleus with low barrier 
and large surface volume ratio (about 70% the nucleons are the 
surface). further emphasize the process our choice nucleus with 
has two effects. (i) The barrier effect combined with high proton 
threshold allows the neutron cross section turn over before the proton cross 
section becomes appreciab'e and thus brings out the difference the (y, 
and (y, processes. (ii) Since would expect even fewer protons 
than usually the case, and thus the presence large number protons 
points the anomaly. From statistical theory would expect that the 
number emitted protons would best equal the number emitted neutrons 
high energy. Our results however show that greater than and 
rising value region where falling with increasing energy. 

The results reported this paper and those (I) silicon and magnesium 
show that the photonuclear process, least these nuclides, cannot ex- 
plained the basis simple evaporation theory, nor does seem that 
minor modifications existing theory will suffice. 


The authors are indebted the National Research Council Canada for 
financial support and for awarding bursary one them (D. Mc). They 
would also like thank Dr. Halpern for transmitting the measurements 
the (y, cross section prior publication. 
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IMAGE INTERFERENCE CALM, NEAR-ISOTHERMAL 


FREDERICK SANDERS 


WITH APPENDIX 


MIRROR” INTERFERENCE THE PRESENCE SLIGHT 
UNIFORM REFRACTION—AN APPROXIMATE THEORY 


STEWART 


ABSTRACT 


Some experiments are described which strong image interference under- 
water sound has been observed and measured. These effects occur frequently 
coastal British Columbia waters, where calm surface conditions are common, 
and have been observed frequencies close kc./sec. The observations are 
compared with theoretical predictions for small sound velocity gradients and 
reasonably good agreement obtained. 


INTRODUCTION 


The phenomenon image interference familiar one the fields optics 
and radio. the former, the classical experiment, wherein 
beam light striking mirror grazing incidence gives rise interference 
fringes, demonstrated nearly all classrooms where elementary physics 
taught. the latter, image interference was not given much attention until the 
use radiation wavelength meters and less became common. case which 
has been much publicized was that the British radar installations (1) 
which guarded the United Kingdom the early years World War II. The 
projection radio energy over relatively flat and highly-reflecting ocean 
resulted very complex lobal structure, and vast efforts were devoted 
methods order eliminate the blind spots produced the 
interference. more recent years radio-astronomers (3) Australia and other 
countries have made good use the phenomenon produce very narrow 
beams for the study galactic noise. 

the case air-borne sound, image interference less frequently observed 
since large areas high reflectivity are not commonly encountered and notes 
constant pitch are rarely employed. When, however, sound projected 
horizontally beneath the surface the ocean, clear that interference may 
expected the frequency constant and the water surface appears 
reasonably flat the sound radiation employed. The existence such phe- 
nomena discussed The Physics Sound the Sea (4). section 
this reference reported that image effect usually marked frequencies 
the order few hundred cycles but that not expected important 
frequencies above 1000 cycles. Young (6) has reported experiments which 
marked interference has been observed frequencies few kilocycles but, 
under the experimental conditions encountered, did not observe any notable 
effects higher frequencies. 


received May 20, 1954. 
Contribution from the Defence Research Board, Ottawa, Ont. 
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Most the work described References and was carried out the deep 
ocean where, even under conditions low zero wind force, good deal 
swell usually exists. the deep enclosed waters the British Columbia 
coast, however, areas flat calm water are often encountered and the 
purpose this paper describe remarkably consistent image interference 
phenomena which have been observed frequencies near the limit the 
audible range. The experiments were carried out Nodales Channel which 
highly-sheltered, over 100 fathoms deep, and through most the year 
characterized nearly isothermal water. The oceanographic features this 
area have been described detail Saur (5). 


THEORETICAL CONSIDERATIONS 


obvious that pulse sound starting point located between the 
surface and the bottom the ocean will have choice several paths 
which can travel receiver also the body the water. The solid lines 
Fig. represent two sound rays which start from point body iso- 


SURFACE 


Fic. Schematic representation sound paths between two underwater points when the 
surface plane reflector. the source, the receiver, and the source image. The under- 
water dotted lines represents the paths under condition slight upward refraction. 


velocity water and proceed point One ray takes the direct path 
while the other arrives via surface reflection The bottom assumed 
too distant too absorbent permit reflection. using the artifice 
virtual image simple matter derive the time delay between 
successive arrivals the same sound pulse via the direct and the surface- 
reflected paths. assume that the separation and large compared 
the depths and obtain for 


where the velocity sound. Recalling that there phase retardation 
radians reflection, clear that there will complete interference 
whenever the range satisfies the relation 

[2] 


where the frequency and integer. Interference will constructive 
when odd, destructive when even. This image interference discussed 
some length Ref. and theoretical sound pressure amplitudes are derived 
terms the range, source and receiver depths, and the coefficient surface 
reflection. 
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While the isovelocity condition valuable introduction the prob- 
lem, does not occur even isothermal water since the velocity sound 
function pressure. water whose velocity varies with depth the sound paths 
will curves whose exact nature will depend the character the velocity 
gradient. The point reflection will displaced indicated the curved 
dotted lines Fig. The problem formulating the law governing interfer- 
ence under these conditions good deal more difficult and has yet not been 
rigorously solved even for the simple case pure linear gradient. Young (6) 
has suggested simplification which makes analytical expression possible 
the case curved ray paths. bases his treatment the argument that 
the difference transit over the two paths that determines the 
relative phase the two signals arriving the receiver. consequence 
should allowable replace the curved paths straight lines and assume 
that the difference travel time along the straight paths, for sound disturb- 
ance whose velocity function depth only, the same would for 
disturbance travelling over the curved paths under the same velocity 
restrictions. This approach yields relatively simple expressions for the time 
difference and, consequently, for the sound intensity distribution. The time 
difference becomes: 


where and are defined by: 


other words, the harmonic mean velocity between the surface and the 
transducer and the harmonic mean velocity between the transducer and 


the receiving hydrophone. The factor given 


which for large compared the product reduces 2ab/r. Equation [3], 
expressed terms order interference becomes: 


[4] 


equation The implications equation [4] will discussed relation 
the experimental observations. 

appendix this paper (pp. 614-619) Stewart has undertaken evalu- 
ate the time difference between the two paths for the special but fairly 
common case where the sound velocity versus depth gradient small but 
constant. considers the times transit over rays which are actually curved 
the natural situation. The only additional limitation imposed that the 
angles between the ray paths and the horizontal small. expression 
for the time difference is: 
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where the gradient divided the surface velocity the harmonic mean 
and the expression expanded, Young’s formula for small constant gradient 
becomes: 


raab 


hence apparent that only for ranges which are very large compared the 
transducer and receiver depths does Stewart’s expression differ materially from 
that Young. 


EXPERIMENTAL PROCEDURES 


September 1952 operation was undertaken obtain accurate measure- 
ments the phenomenon image interference, which was known occur 
very striking form Nodales Channel. During the war years Tully and 
Cameron the Pacific Oceanographic Group, Nanaimo, B.C., had 
noted and made some measurements the effect but security restrictions 
prevailing that time prevented their publication. 

The present experiments were carried out from H.M.C.S. Cedarwood and 
C.N.A.V. Ehkoli, which are the two West Coast research vessels operated for 
the Defence Research Board the Royal Canadian Navy. taut 
between two heavy anchor clumps was first laid Cedarwood that Ehkoli 
could secured fixed position without drifting yawing. Cedarwood was 
the transmitting ship and made closing runs Ehkoli from ranges 2500 
yd. close was deemed safe. Sound frequency 18.0 kc. was pro- 
jected from depth approximately thirteen feet and the sound projector 
trained Ehkoli throughout the run. Sound transmissions were emitted 
rate five two seconds, which 5-knot approach speed amounts 
approximately one pulse per yard. The receiving elements were three hydro- 
phones lowered single steel wire from Ehkoli nominal depths 20, 50, 
and 100 ft. The depth the 100-ft. hydrophone was indicated directly 
repeating depth meter designed Moon Pacific Naval Laboratory. 
Its error, checked still water Esquimalt and again Nodales Channel, 
was less than one foot over the range used. The depth the shallow hydro- 
phone was measured length light silk trout line, drawn taut. This hydro- 
phone was visible all times and, even when the hydrophone was streaming 
badly the tidal current, was reasonably certain that the trout line was 
vertical when depth measurements were made. The depth the 50-ft. hydro- 
phone was obtained plotting the streaming curve the supporting wire 
from the data obtained above. 

The recording equipment consisted three logarithmic amplifiers, four- 
channel oscillograph, and Fairchild oscillograph camera. The fourth channel 
was used record i.e., radio pulse from Cedarwood synchro- 
nized with the sound transmissions. The depth Cedarwood’s sound projector 
was carefully determined both with the ship rest and when was steaming 
knots. From the above experimental arrangements, felt that projector 
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and shallow hydrophone depths were known within in. and the deeper 
hydrophone depths within ft. Range accuracy was +50 ft. and the fre- 
quency was known 100 cycles better. 

Bathythermograph casts were made from Cedarwood, one for each run, the 
location being alternated between the two ends the course. Salinity stations 
were taken three times September 10th and twice September 11th. 
Eleven runs were made the first day and ten the second. The weather 
remained ideal throughout the two days, being overcast with occasional light 


Photograph taken Nodales Channel demonstrate the very calm surface con- 
ditions common these enclosed waters. 


showers and entirely free wind. Fig. photograph taken Nodales 
Channel and representative the very calm surface conditions which 
frequently prevail. 

ANALYSIS RESULTS 


The runs averaged 7000 ft. length that, for each run, over 2000 sound 
pulses were received each hydrophone. desired give the reader some 
idea the regular manner which the sound intensity varies from pulse 
pulse. For this purpose 200 yd. section Run with the 100-ft. hydrophone 
has been plotted Fig. The range this minimum was 1700 ft. but, for 
convenience, transmission number rather than range was used the abscissa. 
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SOUND INTENSITY DECIBELS 


2180 
SOUND TRANSMISSION NUMBER 


observed sound intensity versus the number the transmission (outgoing 
pulse). The intensity greater the right-hand end because the range decreasing. 


With very few exceptions, this type consistency was observed all runs. 
Since the accuracy reading the trace amplitudes the 35-mm. film the 
Fairchild camera was probably not better than one decibel, seen that the 
variability the sound field from one pulse the next was usually only 
slightly worse than the observational error. 

The first step the analysis was study the variation hydrophone depths 
and velocity gradients over the period the runs. Fig. shows the depth 


RUN NUMBER 
= 


HYDROPHONE DEPTH FEET 


SEPTEMBER 10 ——————>| SEPTEMBER 
PACIFIC STANDARD TIME 


Fic. Plot the depth the nominally 100-ft. hydrophone versus time for September 
and 11, 1952. The shaded sections represent the times and durations individual runs. 


the nominally 100-ft. hydrophone plotted against time over the operating 
period. The shaded sections represent individual runs. Around 1300 hr. 
September 10th and again around 1500 hr. September 11th, the streaming 
due tidal current was quite bad. Runs through and through were 
discarded being impossible intercompare. Runs through September 
10th show steady hydrophone depth but examination the bathythermo- 
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graph (BT) traces for this period showed them have variability structure 
the order four five tenths degree between traces taken half-hour 
apart. This was not expected and attributed slightly different 
water drifting through the channel during the water period. The 
period from 0910 through 1320 September 11th, however, showed reason- 
ably satisfactory conditions from the surface 120 ft., demonstrated 
the traces Fig. will noted that all but two traces the water 


DEPTH FEET 


49 49 49. 49 49 49 
50 50 50 50 50 50 50 


WATER TEMPERATURE DEGREES FAHRENHEIT 


Fic. Bathythermograph traces for September 11, 1952. Surface temperatures vary from 
49.2° 49.4° 


was isothermal ft. From 100 ft. there usually slight positive 
gradient* the order one tenth degree Fahrenheit per 100 ft. While 
these gradients are very small, will shown later that they can have signi- 
ficant effects the interference patterns. must noted also that each 
trace represents only the temperature structure along one vertical line the 
2000 yd. water covered the corresponding run. Since was not pos- 
sible obtain more extensive coverage, seems reasonable assume 
that these seven traces represent the probable spread the thermal gradients 
which existed over the sea paths concerned for the period 0910 1320 
September then possible calculate the greatest and smallest 
velocity gradients which occurred over the four-hour period which runs 
through were carried out. 

The salinities observed various depths over this same period were plotted 
against the corresponding temperatures and relation 
parts per thousand per degree Fahrenheit was established. Using Kuwa- 
hara’s tables (2), the maximum and minimum values the velocity gradient 
due temperature, salinity, and pressure were computed for each hydrophone 
depth. From the gradients, the velocity difference equation [4] was 


*In oceanography positive temperature gradient defined mean that the temperature increases 
the depth increases. The same convention adopted with respect salinity and sound velocity. 
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obtained. Table lists the results. Since the absolute velocity enters only 
multiplying factor, the velocities and are listed being the same and equal 
4850 ft. per sec. which was their mean value four significant figures. 


TABLE 
MAXIMUM AND MINIMUM VELOCITY DIFFERENCES 
depth, ft. condition f.p.s. 
Max. 4850 —0.22 
Min. 4850 —0.17 
Max. 4850 —0.93 
Min. 4850 —0.43 
100 Max. 4850 —1.98 
100 Min. 4850 


The actual sound intensities recorded the three hydrophones for runs 
through are shown Figs. and where they are plotted against the 
logarithm the range avoid crowding the short-range end. The great 
variation intensity between peaks and minima apparent nearly all runs. 
Differences db. are common and some cases they amount db. 
When one considers that half-wavelength about in. and that the 
interfering signals have often traversed half-mile water, appears quite 
remarkable that such consistent phase differences and amplitude equalities 
can exist. The individual extrema for the shallow and medium hydrophones are 
readily identifiable from run run. While their ranges are seen wander 
slightly from one run the next, extrema the same order can picked 


SHALLOW HYDROPHONE 


RUN 


RUN 16 


RUN 17 


RUN 18 
RUN 19 
600 1000 200 3000 4000 6000 


RANGE IN FEET 


RELATIVE SOUND INTENSITY —DECIBEL SCALE 


Fic. Plots sound intensity versus log range for the nominally 20-ft. hydrophone. 
Note the rather weak but consistent minimum between 6000 and 7000 ft. This believed 
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MEDIUM HYDROPHONE 


RELATIVE SOUND INTENSITY-DECIBEL SCALE 


RANGE IN FEET 


Fic. Plots sound intensity versus 
log range for the nominally ft. hydro- 
phone. Since the ordinate decibel 
scale some the minima have intensity 
less than one ten-thousandth that the 
adjacent 


RELATIVE SOUND INTENSITY DECIBEL SCALE 


RANGE IN FEET 


Fic. Plots sound intensity versus 
log range for the nominally 100-ft. hydro- 
phone. The orders interference range from 
for the extreme left-hand that 
major minima exist between 4000 and 


7000 ft. 


out casual inspection. Those for the deep hydrophone are more numerous 
and identification little more difficult. The question the order inter- 
ference the various extrema will shortly discussed relation the 
theoretical considerations. 

Before comparing the observed results with theory, appears worth while 
present physical picture the radiation pattern sound projector 
operating beneath calm surface under various conditions refraction. 
the locations the minima only, for simplicity, are displayed 
plot, from equation [2] clear that for zero refraction, they will lie series 
straight lines shown the solid lines Fig. Each line corresponds 
even integral value i.e., path difference even number 
half-wavelengths. Between these lines complete destructive interference 
will lie corresponding group ‘‘lobes’’ sound intensity whose axes would 
also lie straight lines obtained assigning odd values equation 
sound receiver fixed range were lowered through this radiation 
‘pattern, would indicate series sharp minima spaced equal intervals 
depth with strong maxima between them. Alternatively, receiver 
fixed depth were brought slowly from great range, would indicate 
series maxima and minima starting with maximum order which 


would spaced range according the relation 


_ 
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Fic. Loci the minima for 12.9-ft. projector. The solid lines present the isovelocity 
condition for minima orders while the dotted curves refer the minima order 
and under upward refraction due pressure alone. 


now the velocity sound made change with depth, the straight-line 
loci Fig. will become curved. For negative velocity gradient the curves 
will bend downwards and, considering destructive interference only, receiver 
closing the projector from great range will again indicate minima orders 
etc. this case the relation will more complex than that 
equation [2], but the ranges the minima will single-valued and must 
start with order If, however, the gradient positive one, the picture 
becomes great deal more interesting indicated the dotted lines Fig. 
The two curved lines marked and represent the loci minima order 
and for velocity gradient due pressure alone. They were 
obtained inserting equation [4] the values: 12.9 ft., 4850 f.p.s., 
0.0182 

will noted that receiver ft. depth will detect two minima order 
approaches from great range and that the short-range minimum occurs 
noticeably greater range than that predicted for isovelocity water, even 
though the velocity gradient slight. receiver operating ft. will detect 
minima order but should record two order one goes deeper 
receivers the same effect will observed, viz., disappearance the minima 
lower orders and the occurrence two minima for each higher order. keep 
the figure simple, refracted loci have been plotted for the higher orders, but 
essentially the same type upward curving occurs. The loci Fig. were 
obtained the use Young’s expression. For ranges 5000 ft. less and 
receiver depths not over 100 ft., expression (equation [5]) gives curves 
which lie very close those shown. larger ranges Stewart’s curves turn 
considerably more sharply, thus predicting occurrence the second minimum 
shorter range than Young’s expression. will seen later, however, the 
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interference effect long ranges sensitive slight refractive changes 
that successful prediction does not appear too hopeful. 

great deal thought has been given finding simple way comparing 
experimental observations these complex phenomena with theoretical pre- 
dictions. Young’s discussion his observations the major emphasis has 
been comparing the observed and theoretical ‘‘sound i.e., the 
departure the sound intensity from that due simple spherical spreading. 
consequence, sound intensity versus range versus some parameter related 
range has been plotted and discussed. the high frequencies employed 
the Pacific Naval Laboratory experiments, the rate change sound inten- 
sity with range often very great, i.e., maxima and minima occur very close 
together. Any attempt plot observed and computed intensities against 
range the same figure results very confused picture. was accordingly 
decided deal only with ranges and orders interference. Fig. 10, the 


SHALLOW HYDROPHONE 


INTERFERENCE 


MAXIMUM REFRACTION 


ORDER 


MINIMUM REFRACTION 
4 6 8 10 


RANGE THOUSANDS FEET 


Fic. 10. Plot log versus log for the nominally 20-ft. hydrophone. The solid circles 
are mean observed ranges for runs through 19, while the bars represent the observed spread 
range. The dotted curves are from Young’s expression for greatest and smallest depths and 
refraction. With refraction the minima should lie the 45-degree straight lines which 


the curves are asymptotic. 


logarithm order interference has been plotted against the logarithm the 
range. readily seen that equation [2] will yield straight line slope 
degrees when integral values are plotted against the log—log 
scale. The intercepts this 45-degree line will vary with the parameters the 
situation. Equation will yield some type curved line with minimum 

The two dotted curves Fig. represent equation [4] plotted for the 20-ft. 
‘hydrophone with two sets parameters. will recalled that, the dis- 
cussion the oceanographic variables, temperature and salinity, was con- 
cluded that the maximum and minimum refractive conditions should con- 
sidered. Inasmuch there was assignable uncertainty the depths 
and the projector and hydrophone appeared wise plot equation [4] 
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for (i) maximum and with maximum refraction and (ii) 
minimum values and with minimum refraction. Hence the dotted 
curves define the limits within which the range extremum given order 
observed the shallow hydrophone might expected lie. 

Looking back Fig. which gives the observed interference patterns for 
the shallow hydrophone, two very strong minima are observable mean 
ranges 870 and 2000 ft. These are almost certainly the minima 
respective order and less pronounced but quite consistent 
minimum observable about 6000 ft. This assumed the second 
minimum order The solid circles Fig. represent the mean ranges 
the three observed minima while the horizontal bars represent the spread 
the observed ranges. broad way the agreement good, particularly when 
one reflects the small magnitude the path difference required produce 
interference and the relatively great ranges concerned. The predictions 
Young’s theory are quite well substantiated and, for practical purposes, 
would appear that equation [4] useful expression for determining the 
general nature the interference pattern for shallow transmitter and 
receiver. 

Fig. the dotted curves represent the limits within which the minima 
observed the 50-ft. hydrophone might expected fall equation [4] 


MEDIUM HYOROPHONE 


ORDER OF INTERFERENCE 


MINIMUM REFRACTION 
3 4 6 8 10 


RANGE THOUSANDS FEET 


Fic. 11. Plot log versus log for the nominally 50-ft. hydrophone. The reader may 
test other choices order interference displacing the observed points down 
two integral units. other choice asymptotic 45° line with the correct intercept. 


valid expression for the ranges the various extrema under the conditions 
slight refraction observed. The solid circles represent the mean ranges the 
five minima shown the curves Fig. while the horizontal bars represent 
the total spread the observed data. The assignment order interference 
more difficult this case than was for the shallow hydrophone. has been 
assumed that the minimum with mean range 2700 ft. the 
minimum order and that the remainder are orders 10, and 12, 
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respectively. The assignment order the 2700-ft. minimum gives 
unreasonably low log curve and the assignment while 
little more realistic, places the whole series too high. unfortunate that 
there way positively identifying the order interference but, lacking 
such identification, the assumption made above seems reasonable and gives 
quite good fit equation [4]. will noted that, for minimum refraction, the 
second minimum order would expected occur range 10,000 
ft. so. Since the ship runs usually started range 6000 7000 ft., 
data are available this longer range. maximum refraction had existed, 
minima order would expected. 

The most spectacular and, unfortunately, the least tractable set minima 
were those observed the deep hydrophone and displayed Fig. every 
run strong and consistently broad minimum occurred about 3600 ft. This 
appeared most likely the minimum order and the 
solid circles Fig. are plotted this assumption. The open circles represent 


DEEP 


ORDER INTERFERENCE 


MINIMUM REFRACTION 
810 
RANGE THOUSANDS FEET 


Fic. 12. Plot log versus log for the nominally 100-ft. hydrophone. The solid circles 
represent the choice for the 3500-ft. minimum while the open circles are for 
method known for positively identifying the order interference. 


the same data the assumption that 3600 ft. There little choice 
between these assumptions and, consequently, firm conclusion regarding 
the fit equation [4] can drawn. 


DISCUSSION 


The extent which the ranges the observed extrema can compared 
with those predicted theory dependent positive identification the 
_order interference. means for ensuring such identification are known 
exist for the variable-range type experiment just described. the case 
the data for the shallow hydrophone, however, there seems little likelihood 
ambiguity. Referring once again Fig. 10, can seen that any assign- 
ment order other than that chosen would quite unrealistic and appears 
reasonable discuss the implications this figure. 
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While only three points have been plotted, these represent the mean ranges 
the minima observed eight runs. The total spread ranges for each 
minimum quite small and consequently felt that considerable confidence 
can placed these values. Rather less confidence can placed the 
accuracy the temperature and salinity data which determine the velocity 
gradient. There however strong evidence that the lower dotted curve repre- 
sents the very minimum refractive effect possible, since was computed for 
zero temperature and zero salinity gradients. hence represents only the 
contribution pressure, viz., increase velocity 0.0182 ft. per sec. for 
each foot depth. From the oceanographic data taken are confident that 
measurable negative temperature gradient existed over the period runs 
through 19. Even greater confidence felt the absence any negative 
salinity gradient. Negative salinity gradients occur only regions very high 
evaporation rate and British Columbia waters, with extensive fresh water 
being contributed run-off throughout most the year, positive salinity 
gradients are the rule. hence highly probable that the minimum refractive 
condition was greater, rather than less, than that represented the lower 
dotted curve Fig. 10. 

Assuming that the above arguments are valid, interesting note 
that the ranges the observed minima fall consistently rather short those 
predicted theory for even the minimum refractive condition. Fig. 11, 
which displays the data for the 50-ft. hydrophone, the same effect observed. 
Because the higher interference orders involved with the deeper hydrophone, 
there little more possibility wrong assignment order but the choice 
shown seems the most reasonable and again results the observed ranges’ 
falling very close the curve for minimum refraction. 

The first explanation for the low values the observed ranges which occurs 
that Young’s theory may exaggerate the effect refraction. Stewart’s 
approach however essentially rigorous for very small gradients and, 
ranges 3000 ft., his series gives curves almost identical with those 
Young. hence appears that fair amount confidence may placed the 
validity the theoretical expressions the ranges, depths, and gradients 
concerned. possible explanation the discrepancy may found 
natural phenomenon which has far not been discussed, namely the motion 
the water itself. the water were moving the direction sound propaga- 
tion the surface but its velocity decreased reversed with depth, the effect 
would equivalent slight negative sound gradient. difference one 
knot between the surface and 100 ft. would just about equal the normal 
pressure effect. The group runs, through 19, which has been discussed 
was chosen correspond period low current flow, that hydrophone 
streaming would minimum. There is, however, good deal evidence 
the effect that, under stronger tidal current conditions, surface water flow 
and deep water flow are not the same. particular, during periods tidal 
flooding, the deep water has frequently been observed flow quite strongly 
towards the north while the surface water was stationary moving slowly 
southwards. general, the trend Nodales Channel north-to-south 
flow the surface water most tidal periods, with deep south-to-north cur- 
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rents during flood periods. Runs through were carried out during the 
and early ebb periods and, any velocity shear existed, would 
velocity with depth. Since the direction propagation was north-to-south 
all runs, this shear would result equivalent negative gradient. The 
omission quantitative measurements the tidal current profile was ad- 
mittedly most unfortunate. was, however, only very recently that the 
writer and his associates began appreciate the probable significance current 
gradient relation sound refraction. Lacking firm data the tidal current 
velocity profile, all that can said that some downward refraction due 
water movement most probable and that, while this effect would small 
compared normal refraction, would bias the ranges the extrema the 
direction observed. 


CONCLUSIONS 


The conclusions which may drawn from the above experiments are the 
following: 

(i) Under conditions upward refraction, the image interference patterns 
are changed very markedly from those predicted for the isovelocity condition. 
For pressure refraction alone the change great that the simple isovelocity 
formula should never used predict the location the extrema. 

(ii) predicted Young and Stewart, for conditions upward refrac- 
tion, the path difference between the direct and the reflected sound rays 
passes through minimum value. 

(iii) There are two results from this minimum path difference. Certain 
orders interference disappear entirely and the extrema higher order will 
occur two different ranges. limiting condition can occur where the depth 
the receiver and the refraction are such that these two ranges coincide. 

(iv) For relatively shallow projector and receiver depths and ranges 
1000 yd. so, the expressions developed Young and Stewart are capable 
giving the approximate ranges the extrema the velocity gradient small. 
The ranges observed these experiments tend fall short those predicted 
but felt that this discrepancy very probably attributable variation 
tidal current with depth. 

(v) For shallow projectors and receivers and ranges greater than 1000 yd., 
the locations the extrema are sensitive small changes refraction that 
prediction their location not too hopeful. 
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Appendix 


“LLOYD’S MIRROR” INTERFERENCE THE PRESENCE SLIGHT 
UNIFORM REFRACTION—AN APPROXIMATE THEORY 


Introduction 


Young’s theory (6) mirror interference the presence refraction takes 
account the variation the velocity propagation but assumes in- 
correct path that neglects the curvature the rays. While intuitively 
the approximation appears reasonable the curvature small, the limits 
its applicability are not immediately obvious. 

The curvature the path has three effects: 

The actual path length greater than the straight line path. 

The average velocity propagation greater than along the straight 

line path. 

The point reflection displaced the source depth differs from the 

receiver depth. 

exact solution the general case stratified medium appears very 
difficult. However, the simplifying assumption constant velocity 
gradient made, approximate solution may obtained reasonably 
straightforward manner. comparison the results this solution with 
those Young’s theory under the same conditions uniform velocity 
gradient should serve guide the reliability Young’s theory under the 
more general conditions variable velocity gradient. 


Definitions 


shall use Cartesian co-ordinate system: 

increases the range increases from the source 
increases the depth increases from the surface (see Fig. 13). 

velocity the surface. 
and are respectively the sound velocities and 

ay) where constant the velocity gradient uniform. 
angle sound ray with horizontal. 
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Fic. essentially the same schematic drawing Fig. but indicates the nature 
the quantities used develop Stewart’s expression. 


time progression wave front along the ray. 
The surface reflected path designated primed quantities. 


From Snell’s Law: 


where element length along the ray. 


Assumptions 
All angles are small that 


(Note that this implies that large compared with and and the 
approximation breaks down very short ranges.) 

such magnitude that the order of, smaller than, the 
maximum value Therefore and are the order less than the 
maximum value 

assuming that the change cos along any ray negligible terms involving 

That the point reflection may obtained follows: 

direct solution for while theoretically possible since and fully 
determine the position, involves solution rather complicated cubic equa- 
tion, yielding result not readily approximated. shall thus make use the 

device. 


a 
| 
| 
| 
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small reasonable suppose that has value not greatly 
different from the value where which 


ar/(a 
use the first approximation only. 


and 

Thus 


find the time required for wave front travel from the 
direct path, use 


but Snell’s Law any path. Therefore 
cos cos 


where and are the angle and speed respectively the wave front 
assumption 


= A, + ax. 
Substituting equation (iii) and using assumption gives 


maintain the order approximation, sufficient put 
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Hence 


equation (v), the first term the right side represents the time cover 
the horizontal distance; the second the ordinary Lloyd Mirror approximate 
term account for the slant the path; the third term shows the effect 
the variation velocity assuming linear path; and the last term the effect 
the curvature the ray. 


For the surface reflected path: 
Using equation (v) for each section the path, 


Substituting for and equation (ii) gives 


2 
Thus 
(vii) 
2.3 3.5 2 


Each term the right-hand side equation (vii) the same order 
assume that the order a/r (assumption 2). 

The first two terms this expression are the same those given Young 
(6) his result reduced the order approximation used here. The terms 
6)? describe the effect the shift the point reflection, and the 
remaining term gives the effect the curvature the path. 


There are four relevant parameters equation (vii), viz., and 
However, these are not completely independent their effects, and mul- 


(viii) 


may thus plot, for example, curves constant (which are curves 
Three such graphs are shown Figs. 14, 15, and 16. The particular values 
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Fic. This plot analagous Fig. but uses the dimensionless parameters and 
For pressure refraction alone and sound velocity 4850 f.p.s., represents 
projector depth approximately ft. The range corresponding .02 would near 
5300 ft. 


o — 


Fic. 15. Plot versus for For pressure refraction alone and sound 
velocity 4850 f.p.s., this would represent projector depth approximately ft. 


chosen are such that the graphs are representative the effect when 
the order 100 ft. and the same order that produced 
pressure alone sea water. The full lines are the results the whole expression 
equation (viii). 

Except large values the first two terms the expression dominate. 
The effect considering only these two terms shown dotted where differs 
appreciably from the complete expression. The difference between these two 
lines indication the accuracy Young’s approximation. 
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Fic. 16. Plot versus for For pressure refraction alone and sound 
velocity 4850 f.p.s., this would represent projector depth approximately 130 ft. 
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